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Abstract 

We estimate the H i intensity fluctuation power spectrum for a sample of 18 spiral galaxies 
chosen from THINGS. Our analysis spans a large range of length-scales from ~ 300 pc to ~ 
16 kpc across the entire galaxy sample. We find that the power spectrum of each galaxy can be 
well fitted by a power law Pui(U) = A U a , with an index a that varies from galaxy to galaxy. For 
some of the galaxies the scale-invariant power-law power spectrum extends to length-scales that 
are comparable to the size of the galaxy's disk. The distribution of a is strongly peaked with 50% 
of the values in the range a = -1.9 to 1.5, and a mean and standard deviation of -1.3 and 0.5 
respectively. We find no significant correlation between a and the star formation rate, dynamical 
mass, H i mass or velocity dispersion of the galaxies. 

Several earlier studies that have measured the power spectrum within our Galaxy on length- 
scales that are considerably smaller than 500 pc have found a power-law power spectrum with 
a in the range ss -2.8 to -2.5. We propose a picture where we interpret the values in the range 
« -2.8 to -2.5 as arising from three dimensional (3D) turbulence in the Interstellar Medium 
(ISM) on length-scales smaller than the galaxy's scale-height, and we interpret the values in the 
range * -1 .9 to -1 .5 measured in this paper as arising from two-dimensional ISM turbulence in 
the plane of the galaxy's disk. It however still remains a difficulty to explain the small galaxy to 
galaxy variations in the values of a measured here. 
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1. Introduction 



Scale-invariant structures seen in galaxies are believed to be the result of compressible tur- 
bulence. The source of the turbulence could be proto stellar winds, supernova, rotational shear, 
spiral arm shocks, etc. (see lElmegreen and Sc alo 2004 for a review). A var iety of statistical mea- 
sures have been proposed t o characterize the structures seen in the ISM (see Lazarian and Pog osvan 
2000; [Padoan et : alj|200lfc iBrunt and Heverl 120021: lLazarian and Pogosvanll2004t iBurkhart et al. 



20101) . Of these, power spectrum is most popular. In the case of compressible turbulence the 



power spectrum is expected to have a power law shape, and the slope a of the power law con- 
tains information as to the nature of the turbulence. 

In the case of the galacti c neutral hydrogen (H i ) power spectrum analysis was first done by 



Crovisier and Dickevl (1 1983b who found that the slope a of the power is 2.7 for spatial scales 

of ~ 5 pc to ~ 100 pc. Similar results were found by iGreenl (11993b for slightly larger length 
scales, i.e. ~ 200 pc. The very small scales have been probed using absorption studies against 
extended sources. Pow er law slopes in the ra nge ~ -2.5 to —2.8 h ave been found on scales of 
0.01 p c to 3.0 pc bv iDeshpande et all d2000b and iRov et al I j2oTob . Studies of the H i in the 



LMC dStanimirovic et all 1 1999b and SMC dElmegreen et all 1200 lb also show that the intensity 



fluctuations have a power law spe ctrum. 

Recently, Begu m et al. (2006) [henceforth Paper I] presented a visibility -based formalism for 
determining the power spectrum of external galaxies with extremely weak H i emission. Using 

this formalism they have found that the power spectrum of the nearby faint (Mb 10.9) dwarf 

galaxy DDO 210 is a power law, i.e, P(U) — A U a with a 2.6. This is consistent with the 

values of a measured for our Galaxy and the nearby galaxies like LMC and SMC. In a subsequent 
paper (IDutta et all 120081) [henceforth Paper II] we have used the same visibility based formalism 
to measure the H i power spectrum of the spiral galaxy NGC 628. The power spectrum was 
found to be a power law (on scales of 0.8 to 8 kpc), but the slope a was found to be ~ 1.6, 
i.e. much flatter than that determin ed in the earlier studies. The earlier studies all probed much 
smaller length scales (< 500 pc). Dutta et all ((2008) proposed that the difference arises because 
at large scales the measurements corresponds to two dimensional (2D) turbulence in the plane 
of the galaxy's disk whereas the earlier observations were all restricted to length scales smaller 
than the scale height (or thickness) of the galaxy's dis k where three dimen sional (3D) turbulence 



is expected. This picture was subsequently verified dDutta et al.L l2009bh [henceforth Paper III] 
in the galaxy NGC 1058 whose H i power spectrum was found to be well described by a broken 
power law, with slope a - -1.0 + 0.2 at large length scales (1.5 to 10.0 kpc) and another with 
slope a - -2.5 ± 0.6 at small length scales (0.6 to 1.5 kpc). The transition, which occurs at the 
length scale 1 .5 kpc, w as used to est i mate t he scale height of the galaxy's H i disk as 490 ± 90 pc. 
In a more recent study, Block et al. ( 2010b have observed a similar break in the power spectrum 
determined from Spitzer observations of th e LMC and measur ed the line of sight thickness of 
LMC to be in the range 100 - 200 pc. In IDutta etail d2009ab [henceforth paper IV] we have 
extended our study of the H i power spectrum to a sample of five dwarf galaxies. Paper IV also 
presents simulations which substantiat es our interpretatio n of the two different slopes in terms of 
2D and 3D turbulence respectively. In lDutta et al.l (120 10) [henceforth Paper V] we have studied 
the H i power spectrum of a harassed galaxy from the Virgo cluster, where the power spectrum 
slope of the outer part is different from that of the inner part consistent with harassement being 
dominant at the outer parts. 

In this paper we present the power spectrum of H i intensity fluctuations of 18 spiral galax- 
ies drawn from the THINGS sample. THINGS, The HI Nearby Galaxy Survey, is an H i 
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Table 1: Some parameters of the galaxies used for the power spectrum analysis. Columns 1-9 gives 1) Name of the 
galaxy, 2) and 3) Major and Minor axis at a column density of 10" atoms cm~ 2 , 4) Radius corresponding to major axis 
in kpc, 5) Distance to the galaxy, 6) a verage HI inclinatio n an gle, 7) Star Fo r mation Rat e, 8) HI mass, 9) Dynamical 
mass. References are as follows: (1) IWalter etaP <2008t). (21 Ide Blok et alj <2008h. (3) iBottinelli and Gouguenheiml 
<1973|) (4)lHuchtmeier and Witzell fl979l) , (5) iHuchtmeier and Seiradakisl <1985l) , (6) lKamphuis and BriggsN 19921) , (7) 
ISofud tl996h. 
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21-cm e mission survey of 34 nearby galaxies carried out using the NRAO Very Large Array 
(VLA) dWalter et all I2008Q with an aim to investigate the nature of the interstellar medium 
(ISM), galaxy morp hology, star formation and mass distribution across the Hubble sequence. 



de Blok et alj (|2008) present high resolution rotation curves for the galaxies in the THINGS 



sa mple. Star formation rate and efficiency of the THINGS galaxies are also extensively studied 
in Bigi el et alJ (120081) : iLerov et al.1 (120081) . This is as an ideal data set for a comparative study 
of the scale invariant fluctuations observed in the neutral ISM of spiral galaxies. The rest of 
the paper is organized as follows. In Section|2]we discuss the galaxies in our sample and briefly 
mention the power spectrum estimator used in this study. The results are presented and discussed 
in Section[3] Finally we conclude with our main results in the Section|4] 



2. Data and Analysis 

THINGS provides high angular (~6") and velocity (< 5.2 kms 1 ) resolution VLA observa- 
tions of a sample of 34 nearby galaxies spanning a range of Hubble types, from dwarf irregulars 
to massive spirals. In this paper our analysis is restricted to spiral galaxies with minor axis 



We are indebted to Fabian Walter for providing us with the calibrated H 1 data from the THINGS survey. 
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greater than 6'. Table Q] lists different properties of the galaxy sub-sample that we have ana- 
lyzed. The different columns in Table Q] are as follows: column (1) Name of the galaxy, (2) 
and (3)the H i major and minor axis respectively at a column density of 10 19 atoms cm 2 , (4) 
radius corresponding to the major axis in kpc, (5) distance to the galaxy, (6) the average H i 
inclination angle, (7) the star formation rate, (8) the H i mass and (9) the dynamical mass de- 
rived from the rotation curve. The values for the followin g parameters : distan ce to the galaxy, 
star formation rate and the total H i mass are taken from Walter et al. (2008), whereas values 
for the inclination angles are adopted from de Blok et al. ( 2008b . de Blok et al. (2008) presents 
rotation curves for 13 galaxies in our sample derived from the same data as the ones that we 
are using for the current analysis. We use these rotation curves to estimate the dynamical mass 
of these galaxies. For the rest of th e galaxies, the values of t he dy n amical mass are taken from 
Bottinelli and Gouguenheiml ( 1 973h: 1- uchtmeier and Witzell ( 1979 ): Huchtmeier and Seiradakis 
dl985l) : lKamphuis and Briggsl dl992 ): Sofud ( 1 1996b and then rescaled for the adopted distances 
noted in Table Q] Though the galaxy NGC 4826 satisfies our selection criteria, we have not used 
it in our analysis for reasons which are detailed later. 

For our analysis, we started with the calibrated visibility data prior to continuum subtraction. 
In the standard THINGS pipeline, the continuum is subtracted by fitting a linear polynomial to 
each visibility (i.e. AIPS task UVLIN). This method is computationally fast and works very 
well in many cases. However, any extended source, (or even off axis point sources) would have 
visibilities that oscillate with frequency and are not well fit by a straight line. We hence model 
(using clean components) the continuum emission in the field using the line free channels. This 
contribution of the model continuum to the visibilities at all frequencies was then estimated and 
subtracted out using the AIPS task UVSUB. The resulting continuum-subtracted data was used 
for the subsequent analysis. We have identified the spectral channels that contain H i emission 
with a relatively high signal to noise ratio, and used only these channels for the power spectrum 
analysis. 

We have tested the efficacy of our continuum subtraction by comparing the angular power 
spectrum of the channels with H i emission with the the angular power spectrum of the line-free 
channels. The line-free channels give an estimate of the contributions from the residuals after 
continuum subtraction. Our results, shown later in this paper, clearly demonstrate that the H i 
signal is significantly in excess of the residuals from continuum subtraction. The subsequent 
analysis is entirely restricted to the range where the H i signal dominates over the residual con- 
tinuum, and hence we do not expect the contribution from the residual continuum to make a 
significant contribution to the estimated H i power spectrum. 



2.1. The H i power spectrum estimator 

The HI power spectrum can be estimated either from the deconvolved data cubes or the 
calibrated visibilities. The advantage of using the calibrated visibilities is that one is dealing 
directly with the measured quantities whose noise properties are well understood. On the other 
hand most popularly used deconvolutions involve heuristics, and the noise properties of the final 
deconvolved data cubes are not well u nderstood. Previous works t h at used visibility b ased esti- 
mators of the power spectrum include Crovisier and Dickev ( 1983 ): Greenl ( 1993). and LazarianI 
d 1995b . Our method differs from these in the way in which the noise bias is handled. The noise 
bias may be neglected in observations of Galactic HI where the signal to noise ratio is large. 
In observations of external galaxies where the signal to noise ratio is modest, it is important to 
handle the noise bias properly. Paper I and Paper IV together provide a detailed description of 
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the technique that we use to estimate the power spectrum. However, for the convenience of the 
reader, we briefly present the salient features below. 

We model the H i brightness distribution of an external galaxy as 



m = w(0)[i+6im, (i) 

where we assume that the angular extent of the galaxy is sufficiently small so that we may repre- 
sent angular separations as two dimensional vectors 6 in the sky plane. The term W0)I describes 
the global H i distribution of the galaxy ie. the smooth fall off of brightness with angular distance 
away from the center of the galaxy. This can be approximately modeled as an exponential 

W0) = exp( VT2#/0 O ) (2) 

for face-on disk galaxies (see Paper IV) where 6q corresponds to the angular extent of the galaxy's 
H i disk. We conservatively use the angular extent of the semi-minor axis for 6q when the disk 
is inclined to the line of sight. The term 610) represents the fluctuations in the H i brightness. 
These fluctuations are modulated by the galaxy's global H i profile W0) which we can think 
of as a window function. We assume that the fluctuations 51(6) are the outcome of homoge- 
neous and isotropic compressible turbulence, and hence 61(6) is a statistically homogeneous and 
isotropic random field. It is important to note that the assumption regarding the source of the 
fluctuations (ie. turbulence) is not crucial in our interpretation, it could equally well be some 
other mechanism. However, the assumption that fluctuation 61(6) is statistically homogeneous 
and isotropic plays a crucial role in our analysis. The assumption that 610) is statistically homo- 
geneous implies that the two-point correlation function £ defined as 

= «(&)>. (3) 

depends only only on the separation 61 — 62 ie. £(61,62) = £0\ -62)- The assumption of statistical 
isotropy further implies that £ only depends on the magnitude of the angular separation | 6\ — 62 I 
ie. £0i02) = £(\ @i - 6*2 |). The angular brackets in eq. (0) denotes an average over different 
realizations of the random field 61(1, m). 

The angular power spectrum Pm(U) of the H 1 brightness fluctuations is the quantity of 
interest in this paper. The power spectrum Pni(U) is the Fourier transform of the two-point 
correlation function 

Pmtih = J £0)e- a "°- 3 d 2 6. (4) 

where U is a two-dimensional vector that refers to inverse angular scales. The assumptions of 
statistical homogeneity and isotropy imply that £0) and Pui(U) depend only on 6 and U, the 
magnitudes of the respective vectors. 

Now radio-interferometric observations directly measure a quantity related to Pui(U), viz. 
the visibility, 

«V(^) = j 10) e- i2K ° g d 2 6. (5) 

which is the Fourier transform of the brightness distribution. Here the two-dimensional vectors 
U, referred to as baselines, denote the separation d between pairs of antennas measured in units 
of the observing wavelength A (ie. U = d/A). The antenna separations d are projected onto the 
plane normal to direction of observation, i.e a plane parallel to the sky plane. 
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We return to the relation between U and Phi(U) below, but first we discuss the option of 
using eq. <(5j to make a image of the source, and then determine Pui(U) from this image. As is 
well known, it is possible to invert the Fourier relation (eq. [5}, and use the measured visibilities 
to reconstruct the sky brightness distribution 1(0). The visibilities are typically not available 
for all the U values that are required to reconstruct the galaxy's image. Consequently, image 
reconstruction involves complicated non-linear and non-local deconvolution techniques which 
rely on human judgment to a certain extent. Deconvolution is particularly important in our 
context where we are not interested in any particular visible feature in the galaxys image. On 
the contrary, we are intent on quantifying the statistical properties of the random field 51(6). The 
effect of the deconvolution on the 61(9) is difficult to quantify. Further, the noise in the different 
pixels of deconvolved image are correlated - the actual correlation properties are also difficult to 
quantify. For these reasons, it is non trivial to quantify the uncertainties in the power spectrum 
obtained from the deconvolved images. These problems can be avoided if one directly works 
with the visibilities, where the noise properties are relatively straight forward. We describe such 
an estimator below. We note in passing, that although we do not use the image to estimate 
the power spectrum, we have made images of all the galaxies that we analyze. These images 
were only used to estimate the angular extent of the HI disk and to determine the channels with 
significant HI emission. 

The most straight forward visibility based power spectrum estimator dCrovisier and Dickey , 



1983b is the square of the visibilities (] < V(t7) | 2 ). However, in addition to the H i signal from the 
galaxy each visibility also has a noise contribution. The latter introduces a positive noise bias in 
the power spectrum estimator | 'V(U) | 2 . While it is, in principle, possible to separately estimate 
the noise bias and subtract it out, this is extremely difficult in practice. The noise bias is often 
larger than the H i power spectrum, and the uncertainties in calibration and noise estimates make 
it extremely difficult to reliably remove the noise bias. 

Paper I presents a technique to estimate the H i power spectrum, avoiding the problem of 
noise bias. This technique is based on the fact that the visibilities at the different baselines 
within a disk of radius (tt#o) 1 around a baseline U are all correlated and they can all be used to 
estimate the power spectrum Pm(U). We correlate the visibility 'V(U) with all other visibilities 
^(0 + At?) within a disk |A£?| < (nG T l to obtain the estimator P m (U) = CV0) «V*(£? + At?)>. 
The self correlation of a visibility with itself is not included to avoided the noise bias. The 
estimator is further averaged over different directions U and binned to increase the signal to 
noise ratio. The relationship between the estimator 'Ph^U) and the H i power spectrum Phi(^0 
can be shown to be (see paper IV) 

Pm(0) =1 W(U) I 2 ® Pm(U) , (6) 

where W(U) is the Fourier transform of the window function W(G) introduced in eq. ([TJ. For a 
window function centered at G = and with angular extent Go, we expect | W(Cf) | 2 to be sharply 
peaked at U — and to fall off rapidly for U » ftT 1 . Hence, at large baselines, it is adequate to 
approximate | W(U) | 2 by a Dirac delta function in eq. ©, whereby 

Pmitf) = K Pm(U) (7) 

ie. we may ignore the convolution and !Phi(^) gives a direct estimate of the H i power spectrum 
Phi(U). Numerical studies (Paper IV) show that it is possible to identify a baseline U m such that 
eq. (0 holds for all U > U m , and we may directly interpret the estimator !Phi(£?) as the H i 
power spectrum Phi(U) at baselines larger than U m . At baselines smaller than U m , the window 
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function has a non-negligible effect on the slope of the power spectrum estimator. The value of 
U m is proportional to 9 l , and it also depends on the slope of the H i power spectrum (Paper IV). 
In our analysis, we have separately determined U m for each galaxy in our sample and used only 
the range U > U m for the subsequent analysis. It should be noted that the estimator VmiU) has 
both real and imaginary parts, and we use only the real part to estimate the H i power spectrum. 
We expect the imaginary part to be zero if the various assumptions made in our analysis are 
all perfectly valid. The measured estimator usually has a small imaginary component. This 
arises due to the limitation of the various assumptions and also as a consequence of noise. The 
requirement that the real part of the estimator should be considerably larger than the imaginary 
part serves as an useful check of our method. In the subsequent analysis we have also used this 
requirement to constrain the U range where it is possible to make a reliable estimate of the H i 
power spectrum. 

For each galaxy in our sample, we determine if the estimated power spectrum can be modeled 
as a power law Pui(U) = A U a . The largest baseline U max used while doing the power law fit is 
determined by the criteria that the real part of the estimator should be greater than the imaginary 
part, and also greater than the power spectrum estimated from the line-free channels. We first 
make an initial estimate of the smallest baseline U m from a visual inspection of the estimated 
power spectrum. We use the range U m < U < U max to fit a power law through a^- 2 minimization. 
The best fit a obtained by this fitting procedure and the galaxy's angular extent are both used to 
get a revised estimate for U m (see Figure 4 in paper IV). We iterate through these steps a few 
times till it converges, and this gives us a range of baseline (E/ m ;„ to U mnx ) over which we have 
a power law fit to the power spectrum. To test that the window function actually has a very 
small effect, we convolve the best fit power law with | W{U) | 2 assuming an exponential window 
function. The goodness of fit to the data is also estimated by calculating x 1 for the convolved 
power spectrum. We accept the final fit only after ensuring that the effect of the convolution 
can actually be ignored. The method that we have used to determine the best fit power law is 
only sensitive to the slope of the power spectrum and not the amplitude. A detailed model of the 
window function (as opposed to the simple exponential assumed here) is needed to determine 
the amplitude of the H i power spectrum, and we have not attempted this in the present paper. 
Finally, we have also converted the range of baselines (f/ mI „, Umax) to a range of length-scales 
(Rmax,Rmin) using R = D/U where D is the distance to the galaxy. This provides an estimate of 
the range of length-scales where the power-law fit holds. 

Fluctuations in the H i specific intensity could arise from fluctuations in either the density or 
the velocity. The velocity fluctuations, however, ar e not important if the width of the frequency 
slice is larger than the galaxy's velocity dispersion dLazarian and Pogosvanl 120001) . In our anal- 
ysis we have first averaged the visibilities over all the frequency channels that contain significant 
H i emission, and we then use these to estimate the power spectrum. We have checked that the 
resulting channel width 6v is larger than the typical velocity dispersion for all the galaxies in our 
sample. The channel averaging effectively collapses the signal, and we may interpret the mea- 
sured power spectrum as being proportional to the angular power spectrum of the fluctuations in 
the H i column density. 

The lcr error-bars for the estimated power spectrum are a sum, in quadrature, of the con- 
tributions from two different sources of uncertainty. At small U (i.e. large angular scales) the 
uncertainty is dominated by the sample variance which comes from the fact that we have a finite 
and limited number of independent estimates of the true power spectrum. At large U (i.e. small 
angular scales), it is dominated by the system noise in the visibilities. We have used thes e error 
bars in our x 1 minimization procedure. Details of the error estimation can be found in Dutta 
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Table 2: Result of the power spectrum analysis. Column 1 to 7 gives 1) name of the galaxy, 2) width of the channel used 
to estimate the power spectrum, 3) and 4) the range of U value for which the power spectrum is evaluated, 5) and 6) 
correspoindin length scales and 7) the power law index a with 1 — <x error 
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(kA) 


^max 

(kA) 
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NGC 628 


41.6 


1.0 


10.0 


0.8 


7.5 


-1.6 ±0.1 


NGC 925 


41.6 


1.0 


10.0 


0.9 


9.2 


-1.0 ±0.2 


NGC 2403 


83.2 


0.7 


7.0 


0.6 


4.0 


-1.1 ±0.1 


NGC 2841N 


83.2 


1.0 


10.0 


1.4 


14.0 


-1.7 ±0.2 


NGC 284 IS 


83.2 


1.0 


10.0 


1.4 


14.0 


-1.5 ±0.2 


NGC 2903 


83.2 


0.8 


8.0 


1.1 


11.1 


-1.5 ±0.2 


NGC 303 IN 


41.6 


2.0 


10.0 


0.4 


1.8 


-0.7 ±0.1 


NGC 3184 


41.6 


0.7 


7.0 


1.6 


15.8 


-1.3 ±0.2 


NGC 3198 


83.2 


1.6 


10.0 


1.4 


8.6 


-0.4 ± 0.3 


NGC 3521N 


83.2 


1.0 


17.0 


0.6 


10.7 


-1.8 ±0.1 


NGC 3521S 


83.2 


1.0 


17.0 


0.6 


10.7 


-1.6 ±0.2 


NGC 3621 


83.2 


1.0 


12.0 


0.6 


6.6 


-0.8 ±0.2 


NGC 4736 


83.2 


0.6 


10.0 


0.5 


7.8 


-0.3 ±0.2 


NGC 5055 


83.2 


1.0 


10.0 


1.0 


10.0 


-1.6 ±0.1 


NGC 5194 


83.2 


1.0 


8.0 


1.0 


8.0 


-1.7 ±0.2 


NGC 5236 


83.2 


0.6 


6.0 


0.8 


7.5 


-1.9 ±0.2 


NGC 5457 


83.2 


0.6 


12.0 


0.6 


12.3 


-2.2 ±0.1 


NGC 6946 


20.8 


1.5 


10.0 


0.3 


4.0 


-1.6 ±0.1 


NGC 7793 


41.6 


0.6 


6.0 


0.6 


6.5 


-1.7 ±0.2 


IC 2574 


41.6 


1.8 


10.0 


0.4 


3.3 


-1.7 ±0.3 



d201lh . 

As mentioned earlier, the galaxy NGC 4826 meets all the requirements of our selection cri- 
teria. However, this galaxy has a very bright H i core which makes the window function W0) 
complicated, and it is not possible to conclusively determine U m for this galaxy. We hence not 
include this galaxy in our sample. 

3. Result and Discussion 

Figure[TJshows the estimated H i power spectrum for a particular galaxy NGC 628. For this 
galaxy, we show the result over the baseline range 0.5 kA to 30 kA, even though the observational 
data extends to larger baselines (U > 30 kA). The H i signal (Phi(^O) falls at the larger baselines 
which are also sparsely sampled, and the power spectrum estimated at these baselines is too noisy 
to be useful. Further, the residuals from continuum subtraction (Figure [TJ become comparable 
to the H i signal at the large baselines. For channels with H i emission in the U range 1.0 kA 
to 10.0 kA, the real part of the estimator !Phi(^) is larger than the imaginary part as well as the 
power spectrum measured from the line-free channels. We have restrict the power law fit to the 
U range 1.0 kA to 10.0 kA (i.e, t/ m ,„, U max ). For the remaining galaxies we only show results for 
the U range that is useful for determining !Phi(^0- 

Figures|2][3]and|4]show the estimated H i power spectrum for all the galaxies in our sample. 
The results are also summarized in Table|2]where column (2) to (7) gives the velocity width of the 
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collapsed channel, the baseline range used for the fit, the corresponding range of length-scales 
and the best fit power law index a respectively. 

The angular extent of the galaxies NGC 2841, NGC 3031 and NGC 3521 are comparable 
to the telescope's field of view. For each of these galaxies, THINGS contains two distinct sets 
of observations each with a different pointing center referred to as North (N) and South (S) 
respectively. Each data set has been separately analyzed and reported in Table [2] with the suffix 
S or N to distinguish between them. The power spectrum could be estimated only over a limited 
baseline range for NGC 3031S, and hence we do not consider it here. For both NGC 2841 
and NGC 3521, we find that the slope a estimated respectively from the N and the S data are 
consistent with one another. In the subsequent analysis, for these two galaxies we have used the 
values of a corresponding to the N data. 

In our analysis (Table 13, for most of the galaxies we have been able to model the H i power 
spectrum using a power law Pui(U) - AU a across a baseline range which approximately spans 
a decade in U. We find that NGC 303 IN has the smallest dynamical range with (U m i„, U max ) = 
(2, 10)k/l. We recollect that this galaxy has two data sets, and the above results are only from 
the N data, the H i power spectrum from the S data could not be fitted with a power law over 
a significant U range. In addition to this, there are four other galaxies (NGC 3 198, NGC 5194, 
NGC 6946 and IC 2574) that have a dynamical range less than a decade (Umax/U m i„ < 10). The 
power law fit in NGC 5457 covers the largest dynamical range with (t/ m ,„, Umax) = (0-6, 12) LI, 
which corresponds to the range of length-scales 0.6 kpc to 12.3 kpc. The power spectrum is 
found to have a slope a = -2.2 + 0.1 for this galaxy. This galaxy has the steepest slope in our 
sample, and in fact, it is the only galaxy with a < —2. We also find, the smallest length-scale ( 0.3 
kpc) occurs in NGC 6949 where a = -1.6 ± 0.1, and the largest length-scale (15.8 kpc) occurs 
in NGC 3184 where a = -1.3 ± 0.2. Over the range of length-scales probed in this sample, 
there does not seem to be any correlation between the length-scale and the slope a. However, 
it is perhaps more relevant to interpret length-scales in terms to the size of the galaxy. We have 
used the value of r2$ to characterize the intrinsic sizes of the different galaxies in our sampl^l- 
For each galaxy, the length-scales R = D/U have been expressed in units of ^ and R/r25 has 
been shown at the top margins of the plots in Figures [2] G] and [4] Our results clearly shows the 
presence of scale invariant fluctuations with a power law power spectrum that extends to length- 
scales comparable to the radius of the galaxy's optical disk (Rjrzs ~ 1). In fact, we are able to fit 
a power law all the way to R/fys > 2 in NGC 6946. We note that r25 is a measure of the optical 
disk size, the H i disk size is generally somewhat larger. 

Figure [5] shows a histogram of the estimated power law index a. The distribution, we find, 
is sharply peaked in the range a - -1.9 to -1.5. The sample mean and the sample standard 
deviations are —1.3 and 0.5 respectively. It is interesting that 9 out of 18 galaxies in our sample 
have a in the range a = -1.9 to -1.5 irrespective of the finer details of their morphology. 
This indicates the possibility of a common origin for these scale free structures. Note that the 
power law index a for the galaxies NGC 925, NGC 2403, NGC 3031, NGC 3198, NGC 3621, 
NGC 4736 and NGC 5457 are not in the range -1.9 to -1.5. We have visually investigated the 
column density maps of these galaxies to see if the difference in slope is associated with some 
feature in the galaxy's morphology. We find that the galaxy NGC 4736 has a H i ring near its 
center which may cause a complicated window function. However, none of the other galaxies 
have such structures which may cause a steepening or shallowing of the power spectrum. Note 



2 The C2s values have been adopted from the R3 catalog of galaxies 

9 




0.5 1 2 5 10 20 



U [kX] 

Figure 1 : The real (filled triangle) and imaginary (open circle) parts of the estimator evaluated using the channels with 
HI emission is plotted with their lcr error bars for the galaxy NGC 628. The real part for the line-free channels (star) 
is also shown. The best fit power law is shown by the straight line. Only the baseline range where the real part of the 
estimator is larger than its imaginary part and also the real part for the line-free channels is used to estimate the power 
spectrum. 

that R/r25 spans from 0.02 to 0.3 for the galaxy NGC 5457, and we are actually probing length- 
scales which are quite small in comparison to the galaxy's optical disk. The interpretation in 
terms of 3D turbulence for this galaxy is therefore quite well justified. We may expect to observe 
the transition to 2D at the smallest baselines, however large error-bars and the limited spatial 
extent at this end makes it hard to detect this in the present analysis. All the galaxies in our 
sample have spiral arms which can apparently be a source of large scale fluctuations. We have 
verified using simulations that a model galaxy which has smooth emission, apart from large scale 
spiral arms, does not produce a power law power spectrum. Further, the galaxies in our sample 
all have different morphology, and one would expect their power spectra to have different slopes 
if the entire fluctuations in the ISM were just a consequence of the spiral arms. The fact that a 
large fraction of the galaxies in our sample have a similar value of the power law index makes it 
unlikely that the measured power spectrum corresponds to the fluctuations due to the spiral arms. 

The galaxies in our sample have a wide range of inclination angles '/', ranging from 7° for 
NGC 628 to 74° for NGC 2841. In paper IV we have performed simulations to assess the effect 
of the inclination angle on power spectrum estimation. We found that though the length scale 
range over which the power spectrum can be estimated depends on the inclination angles, the 
power law slope remains constant irrespective of the inclination angle. In order to investigate if 
the measured slope is influenced by the galaxy's inclination angle, we have evaluated the linear 
correlation between the power law index a and the inclination angle of the galaxy ( Figure|6]i. The 
lack of correlation (correlation coefficient 0.25) confirms the expectations from the simulations 
that the dispersion in the values of a is not influenced by the inclination. 

It is important to identify the mechanism that generates the large-scale density fluctuations 
observed in the ISM. Here we consider a few possibilities and try to assess if these mechanisms 
can give rise to the a values measured in our sample. 

At small scales, the energy input from star formation processes is a major driving force of 
the ISM turbulence. Star formation process stir up the ISM and the fluctuations hence generated 
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Figure 2: The real (filled triangle) and imaginary (open circle) parts of the estimator evaluated using the channels with 
HI emission is plotted with l<x error bars for the galaxy indicated in each of the individual panels. The real part for the 
line-free channels (star) is also shown. The best fit power law is shown by the straight line. Only the baseline range 
where the real part of the estimator is larger than its imaginary part and also the real part for the line-free channels is used 
to estimate the power spectrum. The top margin of the eachjpanel shows the length scale corresponding to the baseline 
value. The length-scales are shown in units of the rn of the respective galaxy. The rn values are adopted from the R3 
catalogue of galaxies. 
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Figure 3: The real (filled triangle) and imaginary (open circle) parts of the estimator evaluated using the channels with 
HI emission is plotted with 1<t error bars for the galaxy indicated in each of the individual panels. The real part for the 
line-free channels (star) is also shown. The best fit power law is shown by the straight line. Only the baseline range 
where the real part of the estimator is larger than its imaginary part and also the real part for the line-free channels is used 
to estimate the power spectrum. The top margin of the eachj>anel shows the length scale corresponding to the baseline 
value. The length-scales are shown in units of the r25 of the respective galaxy. The r25 values are adopted from the R3 
catalogue of galaxies. 
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Figure 4: The real (filled triangle) and imaginary (open circle) parts of the estimator evaluated using the channels with 
HI emission is plotted with l<x error bars for the galaxy indicated in each of the individual panels. The real part for the 
line-free channels (star) is also shown. The best fit power law is shown by the straight line. Only the baseline range 
where the real part of the estimator is larger than its imaginary part and also the real part for the line-free channels is used 
to estimate the power spectrum. The top margin of the each panel shows the length scale corresponding to the baseline 
value. The length-scales are shown in units of the r 2 s of the respective galaxy. The r 2 s values are adopted from the R3 
catalogue of galaxies. 
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Figure 5: Shaded regions shows the histogram of the distribution of the power law index a. Note that the histogram has 
5 bins for 18 data points. The sample mean ji and sample standard deviation cr are also shown. 
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Figure 6: Scatter plot of the average inclination angle i with power law index a. The 1 
The value of the linear correlation coefficient r is given at the top right corner. 



cr error bars of a is also shown. 
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Figure 7: Scatter plot of the surface density of the star formation rate (Tabled with power law in dex a. Open circles 
show the values from this paper. We also show the same quantities for the five dwarf galaxies from Dut ta et aT] i2009al) 
using open boxes. The l<x error bars of a are also shown. The left panel shows the SFR per unit area of the optical disk 
of the galaxies. The right panel shows similar quantities with the surface density of SFR calculated per unit HI disk area. 
We also evaluate the linear correlation coefficients for each of the two cases using the galaxies in this paper. These values 
are given at the top right corner of the each panel. The x axis of each plot has a unit of kpc~ 2 , while the SFR values are 
considered to be in units of M Q yr~' . 
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Table 3: The values of r^, (cr), cr 2 5 for six galaxies of our sample. These values are adopted from Tambu rro et alj f2009t) . 



Galname 


ri5 




C~25 


a 




(kpc) 


(km s" 1 ) 


(km s _1 ) 




NGC 628 


10.1 


8.0 


7.7 + 0.4 


-1.6 ±0.1 


NGC3184 


11.7 


10.4 


9.7 + 0.9 


-1.3 + 0.1 


NGC 4736 


5.3 


12.0 


12.1 + 1.4 


-0.3 ± 0.2 


NGC 5194 


11.4 


17.7 


17.0 ± 1.4 


-1.7 + 0.2 


NGC 6946 


9.2 


10.1 


8.6 + 0.9 


-1.6 + 0.1 


NGC 7793 


5.7 


11.4 


9.7+ 1.1 


-1.7 + 0.2 



diffuse to larger scales. If the decay process is scale independent, it can give rise to scale free 
structures as seen in our study. In such a scenario, the galaxies with higher global s tar formation 
rate w ould have a shallower power spectrum. For a sample of irregular galaxies, Iwillett et alJ 
(2005) find that on scales of 10 - 400 pc the power law index of the Ha power spectra becomes 
steeper as the star formation rate (SFR) per unit area increases. In paper IV we reported a weak 
correlation between the SFR and the power law index of the H i intensity fluctuations in dwarf 
galaxies. To investigate if star formation is driving the turbulence we see here, we calculate the 
linear correlation coefficient of a with the disk averaged SFR per unit area of the optical disk as 
well as the SFR per unit area of the H i disk. The result, shown in Figure [7] suggests that there 
is no correlation between the SFR and the power law spectral index at the length-scales probed 
in our study. 

ISM turbulence can also be driven by the kinetic energy of the H i gas. In this ca s e, one 
would expect the H i velocity dispersion (cr) to be correlated with a. Tamburro et al. I d2009h 



have estimated (cr) as well as (T25 for 1 1 galaxies in the THINGS sample. We consider here the 
galaxies with inclination angle smaller than 50° to reduce the effect of the contribution from large 
scale velocity gradients in the observed velocity dispersion. These galaxies are listed in Table [3] 
where columns (2) to (5) give (2) r25 in kpc, (3) (cr) in km s _1 , (4) (T25 in km s -1 and (5) a. Note 
that cr 2 5 represents the velocity dispersion at the radius corresponding to r 2 5, within which major 
star formation takes place, whereas (cr) is an average over the entire disk. The maximum length 
scales probed using the power spectrum analysis presented in this paper is comparable to the 
radius r 2 5 for these galaxies, however, the power spectrum estimates are over the entire H 1 disk. 
Figure [9] shows the scatter plot of (T25 and (cr) with a. We find no correlation between a and cr 
for these six galaxies. Strictly speaking, the velocity dispersion is more likely to be correlated 
with the amplitude of the power spectrum rather than its slope. Since the estimator that we have 
used here is not sensitive to the amplitude of the power spectrum, we can not conclusively rule 
out the influence of the velocity dispersion. 

It is possible that the scale-invariant power-law power spectra measured here may be the 
outcome of the gravitational effects of the galaxy's dark matter halo or the self gravity of the 
galaxy's H 1 gas. In such a scenario, it is possible that a is correlated with the galaxy's total 
dynamical mass Mj v or its total H 1 mass Mhi- Figure [8] shows a scatter plot of a with Md x 
(left) and Mhi (right) respectively. The dynamical mass estimates are obtained either from the 



THINGS data dde Blok et a l.. 2008) whenever available, or from earlier references as discussed 



in Section 2. The Mhi values are taken from lWalter et al. (2008). Our analysis yields very low 
linear correlation coefficients suggesting that the gross gravitational effects of the dark matter 
halo or the H 1 disk are not the responsible for the observed power-law power spectra. 
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Figure 8: Scatter plot of the dynamical mass Mj y and the total HI mass Mm (Tabled with power law index a. The 1 —cr 
error bars of a is also shown. The value of the linear correlation coefficient r is given at the top right corner. 
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Figure 9: Scatter plot of (<r) and o"2s (Table[3J with power law index i 
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In summary, the slope a is not correlated with any of the galaxy parameters that we have 
considered. While 50% of the a values are in the narrow range -1.9 to -1.5, the remaining 
galaxies have a values spread over the broad range -2.2 to -0.3. It is unclear if this large spread 
arises from statistical fluctuations in the distribution of a, or if there is an underlying physical 
mechanism that is responsible for this. 



4. Summary and Conclusion 

We have measured the angular power spectrum of 2 1 -cm specific intensity fluctuations for a 
sample of 18 spiral galaxies drown from THINGS. This is the first comprehensive analysis of the 
power spectra of a moderately sized sample of external spiral galaxies. For all the galaxies the 
estimated power spectrum can be well fitted with a single power law across a reasonably large 
dynamical range. The power-law power spectrum, we find, extends to length-scales as large as 
~ 10 kpc indicating the presence of scale-free structures in the ISM on length-scales that are 
comparable to the size of the galaxy. 

In this analysis we have been able to measure the H i power spectrum over a large range of 
length-scales spanning from ~ 300 pc to ~ 16 kpc across the entire galaxy sample. The power 
spectra, we find, are well fit by power laws indicating the presence of scale-invariant fluctuations. 
Fifty percent of the galaxies in our sample have a measured power law index (slope) a in the 
range -1.9 to -1.5 (Figure |5]). Only one galaxy, NGC 5457, has a slope a = -2.2 which is 
steeper than -2. All the other galaxies in our sample h ave slopes a > —2. 

A large number of earlier studies (summarized in Elmegreen and Scalol 2004 ), which have 



probed small length-scales ranging from 10 pc to 200 pc, find a power law power spectrum 
with slope a x -2.8. This is believed to be the outcome of three dimensional (3D) compressible 
turbulence in the ISM. If the same process extends to length-scales larger than the galaxy's scale 
height, we would expect to see two-dimensional turbulence in the plane of the galaxy's disc. 
Dimensional arguments given in Papers II and III lead us to expect the slope of the 2D density 
fluctuations to be a a -1.8. This, within estimated measurement errors, is consistent with the 
slope that we have measured here for most of the galaxies in our sample. This prompts us 
to believe that the small scales and large scale fluctuations may both originate from the same 
physical process, presumably turbulence. 

Energy input from supernova is believed to be the major driving mechanism for turbulence 
at small length-scales. However, it is unlikely that this mechanism would be able to generate 
large scale coherent structures as seen here. Further, the range of length scales that we probe 
here ranges from ~ 300 pc to 16 kpc which has practically no overlap with the ranges of length 
scale at which turbulence is believed to be operational in our Galaxy. 

Our analysis confirms that the power-law index index has no correlation with inclination an- 
gle, SFR, dynamical or total H i mass of the galaxy. At present we do not have any understanding 
of the physical process responsible for the scale-invariant large-scale fluctuations measured here. 
The mass fluctuations in the galaxy's dark matter halo is an interesting possibility that we plan 
to pursue in future. 

Numerical stud i es of the magneto-hydrodynamic (MHD) turbulence (jBeresnvak et~aL , 2005 



Kowal et al. , 2007 : Tofflemire et al. , 2011; Burkhart and Lazarian , 20121) demonstrate that the 



Mach number of the medium is correlated with the density fluctuation power spectrum index. 
Alternatively, the Mach number can also be estimated from the turbulent velocity dispersion of 
the medium for such MHD turbulence. We are presently investigating possibilities to separate the 
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turbulent velocity dispersion in the H i gas from it's thermal counterpart. It may be possible that 
the H i column density as well as the velocity dispersion provides measure of similar quantities 
on the turbulent H i . However, it is not straightforward to carry over the results of the simulations 
performed for MHD turbulence to the turbulence in the H i gas, which may be characteristically 
different. This can be tested in future studies and give rise to a better understanding of the ISM 
turbulence in general. 
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